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^ I Abstract 

Equal time, point to point correlation functions for spatially separated meson 
currents are calculated with respect to a variational construct for the ground 
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t^J- ■ state of QCD. With exact calculations, vector, axial vector and scalar channels 



show qualitative agreement with the phenomenological predictions, where as 
pseudoscalar channel does not. However, the pseudoscalar correlator, when 
approximated by saturating with intermediate one pion states agrees with 



results obtained from spectral density functions parametrised by pion decay 



constant and < tpip > value obtained from chiral perturbation theory. 
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I. INTRODUCTION 

Quantum Chromodynamics (QCD) in the low energy sector is nonperturbative and the 
vacuum structure here is nontrivial 0. The vacuum structure of QCD has been studied 
since quite some time both with quark condensates associated with chiral symmetry breaking 
as well as with gluon condensates |||3[]. An interesting quantity to study with such a 
nontrivial structure of vacuum is the behaviour of current-current correlators illustrating 
different physics involved at different spatial distances. This has recently been emphasized 
in a review by Shuryak and studied through lattice simulations 0. The basic point is 
that the correlators can be used to study the interquark interaction -- its dependence on 
distance. In fact they complement bound state hadron properties in the same way that 
scattering phase shifts provide information about the nucleon-nucleon force complementary 
to that provided by the properties of the deuteron || . 

We have recently considered the structure of QCD vacuum with both quark and gluon 
condensates using a variational ansatz |7j]. We shall use here such an explicit construct of 
QCD vacuum obtained through energy minimisation to evaluate the meson correlators. 

We organise the paper as follows. In section II we recapitulate the results of Ref. . In 
section III we define and calculate meson correlation functions. In section IV we quote the 
results. Section V is devoted to the study of the exceptional case of pseudoscalar correlator. 
Finally we discuss the results in section VI. 

II. QCD VACUUM WITH QUARK AND GLUON CONDENSATES 

We have considered the vacuum structure in QCD using a variational approach with 
both quark and gluon condensates 0. Here we shall very briefly recapitulate the results of 
the same for the sake of completeness. The trial variational ansatz for the QCD vacuum is 

taken as 

\ vac >= U G U F \0 > (2.1) 



obtained through the unitary operators U G and Up for gluons and quarks respectively on 
the perturbative vacuum |0 >. 

For the gluon sector, the unitary operator U G is of the form 



U G = exp (B G t - B G ) 
with the gluon pair creation operator B G given by 



(2.2) 



B, 



G f = ^ [ f$)a a i(kj*a a i(-$d% 



(2.3) 



In the above a a i(k) are the transverse gluon field creation operators satisfying the following 
quantum algebra in Coulumb gauge @ 



b z^xt 



a a i (k),a" j (k) 
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(2.4) 



with a a i(k) annihilating the perturbative vacuum |0 >. Further f(k) is a trial function 
associated with gluon condensates. 

Similarly for the quark sector we have, 



U F = exp (£ F f - B F ) 



(2.5) 



with 



bJ 



h(k)d i(k)(a- k)c\{-k) 



dk 



(2.6) 
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Here h(k) is a trial function associated with quark antiquark condensates. The operators c> 
and c create a quark and antiquark respectively when operating on the perturbative vacuum. 
They satisfy the anticommutation relations 



"It 



(*),<£(#)*]+ = 8r.8*S(k-%) = [c\ r (k)As(k')\ 



(2.7) 



Clearly such a structure for the vacuum eventually reduces to a Bogoliubov transforma- 
tion for the operators. One can then calculate the energy density functional given as 



e = F(h(k),f(k)) 



(2.8) 



The condensate functions f(k) and h(k) are to be determined such that the energy 
density eo is a minimum. Since the functions cannot be determined analytically through 
functional minimisation except for a few simple cases j§|, we choose the alternative approach 
of parameterising the condensate functions as ( with k = \k\), 

smhf(k)=Ae- Bk2/2 (2.9) 

This corresponds to taking a Gaussian distribution for the perturbative gluons in the non- 
perturbative vacuum. Similarly, for the function h(k) describing the quark antiquark con- 
densates we take the ansatz, 

A' 
tan2h(fc)= ( ^ 2fc2 _ 1)1/2 (2.10) 

Further one could relate the quark condensate function to the wave function of pion as a 
quark antiquark bound state || and hence to the decay constant of pion. 

In Ref. 0] the energy density is minimised with respect to the condensate parameters 
subjected to the constraints that the pion decay constant i n and the gluon condensate value 
~ < G a G aiiv > of Shifman Vainshtein and Zhakarov [|K| come out as the experimental 
value of 93 MeV and 0.012 GeV A respectively. 

The results of such a minimisation showed the instability of the perturbative vacuum 
to formation of quark antiquark as well as gluon condensates when the coupling became 
greater than 0.6. Further the charge radius for the pion comes out correctly (R c h — 0.65 
fm) for a s = 1.28. The corresponding values of A and R of Eq. ( |2.10| ) are calculated to be 
A' min ~ 1 and R ~ 0.96 fm. 

With the structure of QCD vacuum thus fixed from pionic properties and SVZ value we 
consider the meson correlators in the next section. 

III. MESON CORRELATION FUNCTIONS 

Consider a generic meson current of the form 



J(x)=r a (x)T a ^(x) (3.1) 

where a; is a four vector; a and /3 are spinor indices; i and j are flavour indices; T is a 4 x 4 

matrix (1,7 5 ,7m or 7m7s) 

Because of the homogeneity of the vacuum we define the conjugate current to the above 
at the origin as, 

J(0) = #(0)1^(0) (3.2) 

with r = 7or^7 

The meson correlation function for the above currents is defined as, 

R(x) =< TJ{x) J(0) > vac (3.3) 

From now on we assume that expectation values are always with respect to the nonper- 
turbative vacuum of our model, hence we drop the subscript vac. 
Hence with Eqs. Q, (O) and Q we have 



R(x) = T aP T' xs < T^ a (x)^(x)iP{(0)4(0) > (3.4) 

This reduces to the identity 

R(x) = T aP V' xs < T%(x)i>{(Q) >< W a (x)4(0) > (3.5) 

The above definition of R(x) is exact since the four point function does not contribute. In 
fact, in the evaluation of Eq. ( |3.4|) we shall have a sum of two terms. The first is equivalent 
to the product of two point functions which is Eq. ( |3.5| ). The second term arises from 
contraction of operators at the same spatial point, related to disconnected diagrams and 
thus can be discarded. 

In Eq. ( |3.5| ) the first term can be identified as the interacting quark propagator 

S(x) =< Tip j (x)ip j (0) > 



It can be shown using the CPT invariance of the vacuum [0 that the second term is given 



as 



< T^(x)#(0) > = -75^(x) 7 5 
= -S(-x) 



(3.6) 



Hence the correlation function of Eq. (|3.4j) becomes 



R(x) = -Tr S{x)T'S{-x)T 



Similarly the correlator for massless noninteracting quarks can be given as 



(3.7) 



Rn(x) = —Tr 



S Q (x)r' S (-x)T 



(3.8) 



Our task is now to evaluate the expression (|3.71 ) with the ansatz for QCD vacuum as 
given in Eq. (|2.10|) . Further we shall be interested in evaluating the equal time point to 
point correlation functions. With this in mind we first calculate the equal time interacting 



Feynman propagator in this nonperturbative vacuum given as [T2|] 

S a0 (x) -- 



' %(*),$ (0) 



(3.9) 



Remembering that the expectation value is with respect to the \vac > as given in Eq. 
fcj.lp, the above propagator reduces to 

1 1 



S{x) 



e lk - x dk 



sin2h(k) - (7 • k) cos2h(k) 



2 (2vr) 3 
In evaluating the above we have used the expectation values 



(3.10) 



< 



C(f) t ^(y)>=(2vr)- 3 ^7(A_(A : )) 



'/3a 



e -ik.{x-y) d ^ 



< ^WM ] >= (2vr)- 3 ^/ (K + {k)) a / HS '^dk 



(3.11a) 



(3.11b) 



where, 



A ± (jfc) = -(l±7°sin2/i(fc)±(a-A;) co$2h(k)) 



(3.12) 



We next use the condensate function as given in Eq. ( |2.10| ) to evaluate Eq. ( |3.10| ). Here 
we shall take the parameters A' = 1 and R = 0.96 fm which give the correct pionic properties 
0. Then the equal time interacting quark propagator Eq. ( |3.10|) reduces to, with x = \x\, 



S(f) " 2^ x± + (2n)^2R3 e ~ X2/(2R2) (2^ x* 1{x) (3 ' 13) 



where 



,>-' 1.2 



/*°° / sin kx \ ke 

I(x) = J o [™kx-—) i + {i _ e _ R2k2y/2 dk (3.14) 

We may further note that when h(k) = i.e the condensates vanish, we recover the free 
massless propagator 

S ( x ) = -l-±- [dkj-ke^ s (3.15a) 



2 (2vr) 3 „ 
i 7 • x 

~2n 2 x 4 



(3.15b) 



Having obtained the propagators, we can calculate the correlation function, Eq. ( |3.7| ) for 
a generic current of the form as in Eqs. ( |3.1| ) and ( ^.2| ). For convenience, we will consider 
the ratio of the physical correlation function to that of massless noninteracting quarks 



*(*) (, , Ijr/.^ , T*\-^/rf ^^ 



r'r 



l + V/(x) +^e"^ 2 — — T i-i— (3.16) 



i?o(x) V 2 v V 8i? 6 i ! iJ'Tr[fr'7JT] 

which is then evaluated in different channels with the corresponding Dirac structure for the 
currents. 

IV. RESULTS 

We have studied the above ratio of correlators for four channels. In each channel we 
associate the current with a physical meson having quantum numbers identical to that of 
the current. The results are shown in Table | and in Fig |l[ We may notice some general 
features and relationships among the correlators. The pseudoscalar correlator is always 
greater than the scalar correlator and vector correlator is greater than the axial vector 
correlator. We may emphasize here that these relations are rather general in the sense that 
they do not depend on the explicit form of the condensate function and arise due to the 
different Dirac structure of the currents which is reflected in the generic expression for the 



correlation functions as in Eq. ( |3.16|) . The behaviour of each channel is consistent with that 



predicted by phenomenology except in the pseudoscalar case where the ratio does not go as 
high as expected from phenomenology. We examine this in the next section. 

V. PSEUDOSCALAR CHANNEL 



The explicit evaluation of the pseudoscalar correlator gives, using Eq. ( |3.16|) 



R(x) , I -, , - n ./'" _,.j ,,2 



1 



1 + -x 2 I(x) 



2 „r -.6 



e~ x /R (5 1) 

8R* e [bA) 



Ro{x) 

which may also be read off from column 4 of Table |. This is plotted as a function of x 
in (Fig. 0). As may be seen from (Fig. [I]) this ratio has a maximum of ~ 1.2 at x ~ 1.3 
fm . Phenomenologically || the peak is at ~ 100 at x ~ 0.5. In order to compare our 
results with other calculations we evaluate the same correlator approximately by saturating 
intermediate states with one pion states. 

Using translational invariance the correlator may be written as 

R(x) = < TJ p (x)J p (0) > (5.2) 

= \ (< ^ P (0) J p (0) > e# M + < J P (0)J P (0) > e~ 1 ^) (5.3) 

Also using the fact that for the pseudoscalar current J p = u^^d and J p = —J p we have 

R(x) = - < J p (0)J p (0) > (e tpS + e~ ipS ) (5.4) 

We now insert a complete set of intermediate states between the two currents but retain 
only the one pion state in the sum for the four point function. Thus, 

R{x) = - f <vac\ J p (0) | n a (p) >< n a (p) | J p (0) | vac > {e ip ^ + e~ ip - s ) dp (5.5) 

We may evaluate the above matrix element using the definition of the pion decay constant 



given as 
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< vac | Jt(x) | *-(p) >= (27r) ^ o)1/2 e^ (5.6) 

8 



where Jt* a = [ip , y tJ- 'j 5 T a ilj] is the axial current. 



It can be shown [[14] that the divergence of the axial current gives the pseudoscalar 
current 

d M [V>y*7 5 T>] = 2 i m q [ipj 5 r a ip] (5.7) 



where m q is the current quark mass. Thus taking divergence of both sides of Eq. (" 15760 and 
using Eq. (|5.7|) we get, 

— f TYl 

2 m q < vac \ iJ~(x) \ n a (p) >= ^^y/^ (5-8) 

where we have used p 2 = m 2 . 

In an earlier paper || within our vacuum model and using the fact that pion is an 
approximate Goldstone mode it was demonstrated that saturating with pion states, gives 
the familiar current algebra result 

Tfl — 

m\ = — -f < VV> > (5-9) 

J ■K 

With this result we eliminate quark mass m q in Eq. ( |5.8|) in favour of the quark condensate 
to get the relation 

< vac | J»(g) | .«(p) >= iM ,; M ,;f ^ (5.10) 

The expression for the pseudoscalar correlator now becomes 

The above integral can be evaluated using the standard integral [ffH] 

r P (p 2 + P 2 y- 1 / 2 sin(ap)dp = ^= C^\ cos(utt)T(u + \)K v+l {aP) 

JO \/7T V Oi I 2 



for a > 0, Ref3 > and in the limit v — > 0. We then finally get for the correlator (using 
saturation of pion states) 

flW = is? (-fc-j — ; — (512) 

9 



The correlator for free massless quarks as calculated in the earlier section for pseudoscalar 

is 

Ro(x) = ^ (5.13) 

Hence the ratio is 

R(t\ TT 1 I <^ ?7ll/l >>\ 

xrn^K^m^x) (5-14) 



R{x) _ n 2 ( < ipip >\ 



Ro(x) 16 V A 
We have plotted in Fig. |2| this ratio for our value of < tpip > and that used by Shuryak 

PJPtH . Note that our value of < ipip > is an output of the variational calculation consistent 



with low energy hadronic properties @. We thus observe that the approximate calculation 
of the pseudoscalar correlator due to saturation with one pion states (Fig. @(a)) yields 
higher values (~ 15 times more) as compared to the exact calculations (Fig. p]). Thus 
the fermionic condensate model for QCD vacuum |7|] does not give as high values for the 
pseudoscalar correlator as required by phenomenological results. The value we have used for 
< ipip >~ (190 MeV) 3 is smaller than Shuryak's value of (307.4 MeV) 3 |16[ which appears 
in the parameterisation of the physical spectral density through the coupling constant JT7 



With his value of < ipip > the ratio R(x) / Rq(x) is shown in (Fig. ^|(b)) which agrees with 
phenomenology. 

VI. SUMMARY AND DISCUSSIONS 

We have evaluated the mesonic correlators in this paper using a variational construct 
for the QCD vacuum. Except for the pseudoscalar channel the results show qualitative 
agreement with phenomenological results 0. Following Shuryak |5|,[T7J1, we also see that using 
current algebra approach the pseudoscalar correlator rises sharply with spatial separation. 
Let us recall that the current algebra result also follows from the approximation of saturating 
by one pion states in the normalisation of the pion state 0. 

It might appear that by suitably changing the value of < ipip > in the exact calculation 
one might be able to reproduce all the phenomenological results. Actually we find that it is 

10 



not so. In fact, it adversely affects the correlators in the other channel which can be seen in 
the exact expressions given in column 4 of Table |I[ 

In view of these findings, it is not clear whether saturation of intermediate states by 
one pion states only in the evaluation of the correlator, is sufficiently well justified. We 
therefore think that a unified treatment of correlation functions in all the channels is still 
not available. 
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TABLES 



TABLE I. Meson currents and correlation functions 



CHANNEL CURRENT 



PARTICLE 



(J P ,MASS in MeV) 



CORRELATOR 

R{x) 




Pseudoscalar 



JP = u l5 d 7r°(0-,135) 



1 + \x 2 I{x) 



I 7T X 6 -X 2 /R 2 



Scalar 



J s =ud 



none(0 + ) 



1 + \x 2 I(x) 



7T X 6 p -X 2 /R 2 



Vector 



J„ = u lf ,d p ± (l',770) 



1 + \x 2 I(x) 



, 7T X 6 f> -X 2 /R 2 



AW 



Axial 



4 = uy^5d ^(l+.llOO) 



1 + \x 2 I{x) 



n x 6 e -x 2 /R 2 



i The integral I(x) is defined in Eq. ( |3.14| ) 
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FIG. 1. The ratio of the meson correlation functions in QCD vacuum to the correlation func- 
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tions for noninteracting massless quarks, — — , Plotted vs. distance x (in fm) 
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Fig. 2a 



Fig. 2b 
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FIG. 2. The pseudoscalar correlator plotted for our value of < ipip >= (190 MeV) 3 in (a) and 
that of Shuryak < ijjip >= (307.4 MeV) 3 in (b) 
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